Steps in the exact time-dependent potential energy surface 
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We study the exact Time-Dependent Potential Energy Surface (TDPES) in the presence of strong 
non-adiabatic coupling between the electronic and nuclear motion. The concept of the TDPES 
emerges from the exact factorization of the full electron-nuclear wave-function [A. Abedi, N. T. 
Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)]. Employing a ID model-system, 
we show that the TDPES exhibits a dynamical step that bridges between piecewise adiabatic shapes. 
We analytically investigate the position of the steps and the nature of the switching between the adi- 
abatic pieces of the TDPES. 
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The description of coupled electron-nuclear motion 
is one of the biggest challenges in condensed-matter 
physics and theoretical chemistry. Fundamental to 
our understanding is the adiabatic separation of elec- 
tronic and nuclear motion embodied in the Born- 
Oppenheimer (BO) approximation. It allows one to vi- 
sualize -approximately- a molecule as a set of nuclei 
moving on a single Potential Energy Surface (PES) gen- 
erated by the electrons in a specific electronic eigenstate. 
The BO approximation breaks down when two or more 
BOPES come close or cross. Some of the most fascinating 
and most challenging molecular processes occur in the 
regime where the BO approximation is not valid, e.g. ul- 
trafast nuclear motion through conical intersections HI , 
radiationless relaxation of excited electronic states [2], 
intra- and inter-molecular electron and proton trans- 
fer [3], to name a few. The standard way of studying and 
interpreting these, so-called, "non-adiabatic" processes 
is to expand the full molecular wave function in terms 
of the BO electronic states. Within this expansion, non- 
adiabatic processes can be viewed as a nuclear wave 
packet with contributions on several BOPESs, coupled 
through the non-adiabatic coupling (NAC) terms which 
in turn induce transitions between the BOPESs. While 
this provides a formally exact description one may nev- 
ertheless ask: Is it also possible to study the molecular 
process using a single PES? This question is particularly 
relevant if one thinks of a classical or semi-classical treat- 
ment of the nuclei where a well-defined single classical 
force would be highly desirable. 

In a recent Letter, we have introduced an exact time- 
dependent potential energy surface (TDPES) that, to- 
gether with an exact time-dependent vector potential 
govern the nuclear motion. These concepts emerge from 
a novel way to approach the coupled electron-nuclear 
dynamics via an exact factorization of the electron- 
nuclear wave function [4J. Features of the exact TDPES 
were studied in the presence of strong laser fields QH5). 
In the present Letter we investigate the generic features 
of the exact TDPES without external laser but in the pres- 
ence of strong non-adiabatic couplings. A major result 



will be that the exact TDPES exhibits nearly discontinu- 
ous steps connecting different static BOPES, reminiscent 
of Tully's surface hopping [6] in the classical limit. 

In H we have proved that the exact solution of the 
time-dependent Schrodinger equation (TDSE), 



H^(r,K,t) = id t *(r,R,i), 



(1) 



of the complete system of interacting electrons and nu- 
clei can be written as a single product (unlike the BO 
expansion), ^(r, R, t) = $R(r, t), of the nuclear 

wave-function, x(R, t), and the electronic conditional 
wave-function, $R.(r, i), that satisfies the partial nor- 
malization condition (PNC), / dr|$R,(r, i)| 2 = 1. In the 
absence of time-dependent external fields, the system is 
described by the Hamiltonian H, 



H = H BO (r,R) + T n (R), 



(2) 



that contains the traditional BO electronic Hamiltonian, 
Hbo fe R) = f e (r) + W ee (r) + V en (r, R) + W nn (R), and 
the nuclear kinetic energy, T n (R) . Throughout this pa- 
per we use atomic units (unless stated otherwise) and 
the electronic and nuclear coordinates are collectively 
denoted by r and R, respectively. 

The exact electronic wave-function satisfies the equa- 
tion 



H BO (±m + UZ up - e(Ri)J%M) = t %fei), 

(3) 

where the electron-nuclear coupling operator 

is UZ up E£i(HV„ A,(R,i)) 2 /2 + 

(-iV„x/x + A„(S,t)) (-<v„ - A„(R<)) )/M v . 
The time-evolution of the nuclear wave-function is 
governed by the Schrodinger equation: 



N„ 



(-iV„ + A^Ri)) 2 
2A/„ 



+ e(R,t))x(R,t)=i9 tX (R,t). 

(4) 



These equations lead to rigorous definitions of the TD- 
PES and the time-dependent vector potential 



e(R,t) = e gl (R,t) + e gd {R 1 t) 
A„(R,t) = ($ R (t) -iV„$ s (t) 



— / r 



(5) 
(6) 



The TDPES consists of two parts: e g i (R, t), defined as 

e gi (R,t) = (^(t) H BO (g,U)+U^ up ((^R,t) $r(*) 



— / r 



(7) 

is form-invariant under the gauge-transformation 

%fet) -> - exp^(E,t))%fci), 

x(E,t) -> x(E,*) = exp(-W(E,t))x(S.*)/ whereas 
e g d(R, t), defined as 



e sd (R,f) = ($ R (t) 



*r(*) 



(8) 



is the part that depends on the choice of the gauge. Here, 
(..|..|..) r denotes an inner product over the electronic 
variables only. 

Why is this representation of the correlated electron- 
nuclear many-body problem exciting? The wave- 
function x(R, t) that satisfies the exact nuclear equa- 
tion of motion j4| leads to an TV-body density T(R, t) = 
lx(B;^)| 2 an d an A^-body current density J„(R, t) = 
Im(x*V v x) + r(R, t)A u which reproduce the true nu- 
clear iV-body density and current density obtained from 
the full wave-function ^(r, R, t) 0. In this sense, 
x(R, t), can be viewed as the proper nuclear wave- 
function. The time evolution of x(R, t), on the other 
hand, is completely determined by the TDPES, e(R, t), 
and the vector potential, A„(R, t). Moreover, these po- 
tentials are unique up to within a gauge transformation. 
This uniqueness is straightforwardly proven by follow- 
ing the steps of the current-density version [6| of the 
Runge-Gross theorem \ 7\. In other words, if one wants a 
TDSE Q whose solution x(E> t) yields the true nuclear 
A^-body density and current density, then the potentials 
appearing in this TDSE are (up to within a gauge trans- 
formation) uniquely given by Eqs. l(5][8}; there is no other 
choice. This also implies, that the gradient of this exact 
TDPES is the only correct force on the nuclei in the classi- 
cal limit (plus terms arising from the vector potential, if 
those cannot be gauged away). The goal of this Letter 
is to find out how this exact TDPES looks like when one 
has strong non-adiabatic couplings in the traditional ex- 
pansion in BO states. One major result will be that the 
exact TDPES shows a nearly discontinuous step when- 
ever the nuclear wave-packet splits in the vicinity of an 
avoided crossing of the BOPES. 

To study the exact TDPES we first of all need a prob- 
lem that is simple enough to allow for a numerically ex- 
act solution and that nevertheless exhibits the character- 
istic features associated with strong non-adiabatic cou- 
plings, such as the splitting of the nuclear wave packet. 
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FIG. 1. Left: The first two BOPESs (indicated in the figure) 
together with the 3rd BOPES (black dashed-line) and the ini- 
tial nuclear wave-function (black solid-line) . Right: Adiabatic 
electronic conditional densities as indicated in the figures 



For this purpose we employ the model of Shin and 
Metiu |8|. It consists of three ions and a single electron. 
Two ions are fixed at a distance of L = 19.0 ao, the third 
ion and the electron are free to move in one dimension 
along the line joining the two fixed ions. The Hamilto- 
nian of this system reads 

1 d 2 
2M0R 2 



1 d 2 



1 



1 



If- -HI 



R\ 



erf 



(9) 



Here, the symbols r and R are replaced by r and R, 
the coordinates of the electron and the movable nu- 
cleus measured from the center of the two fixed ions. 
M = 1836 a.u. and we choose Rf = 5.0 ao, Ri = 3.1 cto 
and R r = 4.0 a a such that the first BOPES, e^L is 



(2) 

strongly coupled to the second BOPES, e BO , around 
the avoided crossing at R ac = —1.90 a and there is 
a weak coupling to the rest of the surfaces. The first 
three BOPES are shown in Fig. [T] (left panel), together 

with the BO conditional electronic densities ^^(Ol 2 
(right panels). As expected, \(jr^ (r)| 2 and 
r ) | 2 exhibit abrupt changes, along the i?-axis, at the 
position of the avoided crossing, R ac : \cj)^' (r)| 2 switches 
from being localized around the fixed ion on the left 
(r = —9.5 ao), to be localized around the one on the 
right (r = 9.5 a ); I^^MI 2 on tne other hand, presents 
a single-peak structure for R < R ac and a double-peak 
structure for R > R ac . 

We suppose that the system is initially excited to 
S 2 ) 

-BO 



and |<4 2) (r)| 2 



and the initial nuclear wave-function is a wave- 



packet with the width a = 1/ a/2785, centered at R = 
—4.0 ao (see Fig. [lj black solid-line), i.e., the initial 

full wave-function is * Q (r,iZ) = Ae-^-^/^^ir) 
with A being a normalization constant. Starting with 
^o ( r j R) as initial state, we propagate the TDSE, nu- 
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FIG. 2. First panel (top): The gauge independent part of the 
TDPES (black solid-line) plotted at four different times (in- 
dicated), (red dashed-line) and eg Q (green dashed-line). 
Second panel (from the top): the gauge dependent part of the 
TDPES is plotted at the same times. Third panel (from the 
top): the exact nuclear density (black dashed-line) is shown 
together with \Fi(R, t)\ 2 (red solid-line) and \ F 2 (R, t)\ 2 (green 
solid-line). Lowest panel: the exact time-dependent electronic 
conditional density, \^r{t, t)\ 2 , is plotted. The color range is 
the same as Fig.^ 

merically exactly, to obtain the full molecular wave- 
function \I>(r, R, t) and from it we calculate, as discussed 
in Ref . [5J, the TDPES in the gauge where the vector po- 
tential is zero. Hence, the TDPES is the only potential 
acting on the nuclear sub-system. In the upper panel 
of Fig. |2j the gauge-invariant part of the TDPES y\, 
e g i, is plotted (black solid-line) at four different times, 
along with the two lowest BOPESs, (red dashed- 

(2) 

line) and e BO (green dashed-line). In the second panel 
(from the top), the gauge-dependent part of the TD- 
PES l[8}, e g d, is plotted at the same times. In the third 
panel (from the top), the exact nuclear density (black 
dashed-line), \x(R,t)\ 2 , is shown together with the ab- 
solute value squared of the projection of the full wave- 
function on the first and second BO electronic states, 
i.e., \Fi(R, t)\ 2 = | / dr(/) { ^*(r)^(r,R,t)\ 2 (red solid-line) 

and \F 2 (R,t)\ 2 = \J drc$* (r)tf(r, R,t)\ 2 (green solid- 
line). In the lowest panel, |$7?(r, i)| 2 is presented. 

At the initial time (t = 0), due to the choice of the ini- 

(2) 

tial state, the TDPES coincides with e B ' . Since *f?o(r, R) 
is not an eigenstate of the Hamiltonian d9), it evolves in 
time. At t = 9.0 fs, e gi coincides with for R < R ac , 
goes smoothly through the avoided crossing region and 
follows e B Q for R > R ac , resembling the diabatic PES of 
state 2 in Ref. |8 1, in which the electron interacts with the 
fixed ion on the right (r = 9.5 cto) and with the moving 



ion, but not with the fixed ion on the left (r = —9.5 ao). 
As e g d is constant in this region (Fig. [2|| 7 the TDPES is 
identical with e g i [9 1. The nuclear wave-packet is driven 
by the TDPES to spread towards the avoided crossing 
of two BOPESs, where a significant non-adiabatic tran- 
sition happens and the exact nuclear density splits. Al- 
ready at this moment, a slight transition of the nuclear 
wave-packet to the lower surface is visible around the 
avoided crossing. At later times, e.g., t = 16.22 fs, t = 
26.24 fs and t — 57.35 fs., far from the avoided cross- 
ing, t g i contains steps that connect its different pieces 
that are on top of different BOPESs in different slices of 
it-space. In the region around R ac , it follows the dia- 
batic surface that passes smoothly through the avoided 
crossing. On the other hand, e g d is piecewise constant 
and presents similar steps as e g i . Therefore, the TDPES, 
£gi + tgd, preserves the features mentioned before, i.e., 
(i) far from the avoided crossing, it presents steps that 
connect the regions in i?-space in which the TDPES has 
the shape of one BOPES to the regions in which it has 
the shape of the other BOPES; (ii) around the avoided 
crossing, it follows the diabatic surface that smoothly 
connects one BOPES to the other. 

The exact TDPES represented in Fig.|2]can be viewed 
from a different perspective. The nuclear wave-packet 
from a semi-classical point of view can be represented as 
an ensemble of classical trajectories, along which point- 
particles evolve under the action of a classical force 
which is the gradient of e ff j. According to our obser- 
vations, on different sides of a step such a force is cal- 
culated from different BOPESs. This is reminiscent of 
the jumping between the adiabatic surfaces in algorithms 
such as Tully's surface hopping IflOl ITT1 . However, 
while Tully surface hopping is a stochastic algorithm, 
the jumps in the exact TDPES correspond to an exact 
solution of the TDSE. When the time-dependent vector 
potential can not be set to zero, a gauge can be chosen 
in which e g d is zero and a time-dependent vector po- 
tential together with t g i specifies the classical force that 
the nuclei experience in different slices of i?-space. In- 
vestigating e g i together with the time-dependent vector 
potential for a wide range of situations may help to im- 
prove the existing semi-classical procedures to simulate 
non-adiabatic nuclear dynamics. 

The exact time-dependent electronic conditional den- 
sity, shown in the lower panels of Fig. [2] at different 
times, behaves similarly to the TDPES: (i) it smoothly 
connects a \cj)^ (r) | 2 -like structure, by crossing R ac , with 
a |</>^(r)| 2 -like structure, or vice versa, presenting a di- 
abatic behavior, e.g. at t = 9.0 fs; (ii) it displays abrupt 
changes, between regions that piecewise match different 
adiabatic conditional densities. 

In order to analyze the behavior of the TDPES, we 
rewrite it by expanding the exact electronic conditional 
wave-function in terms of the adiabatic electronic states 



3 



151 . Due to the choice of the parameter in the Hamil- 
tonian, we only need to include the first two BO states, 
then 



$ R (r, t) = d (R, t)ct><£ (r) + C 2 (R, t)<t>$ (r). 



(10) 



We expand the full electron-nuclear wave-function in 
the same basis, 

*(r,JJ,t) = F 1 {R,t)^\r) + F 2 (R,t)<l>f{r), (11) 

where the expansion coefficients, F k 's and C k 's, are re- 
lated as 



C k {R,t) 



F k (R,t) = e- ie ^F k (R,t) 
x{R,t) ^^F 1 (R,t)\ 2 + \F 2 (R,t)\ 2 



(12) 



Here, 8 is the phase of the exact nuclear wave-function 
and we have used the relation \x{R, t)\ 2 — \F\(R, t)\ 2 + 
\F 2 (R,t)\ 2 , determined by J dr\^(r, R)\ 2 using Eq. JTT 
and the orthonormality of the adiabatic states. By us 
ing Eqs. |5) and (lO) , we rewrite e g i(R, t) and e g d(R, t) in 
terms of e£ (R) and C k {R, t) (k = 1, 2) 

e gi (R,t)= \C k (R,t)\ 2 eP (R) (13) 

fe=l,2 

e gd (R,t)= \C k (R,t)\ 2 %(R,t), (14) 

fe=l,2 

where 71 and 72 are the phases of Ci, C 2 . In Eq. | [T3) , 
all terms of C(Af _1 ) have been neglected and it only 
contains BOPESs which are the leading terms respon- 
sible for the shape of e g i(R,t), especially far from the 
avoided crossing where the NACs are small. The gauge- 
dependent term is written in terms of the time deriva- 
tive of the phases, 71 and j 2 . |Ci| 2 an d |C 2 | 2 vary 
between and 1 and |Ci| 2 + \C 2 \ 2 = 1 by virtue of 
the PNC. Therefore, as Eq. (13) suggests, in the region 
where e gi (R,t) coincides with eg (R), the correspond- 
ing expansion coefficient \C\\ 2 is close to one while 
\C 2 \ 2 is close zero and vice versa. We have observed 



(Eq. 1 12 1) that at R , the cross-over of and \F 2 \ where 

|Fi(i?o,t)| = 1*2(^0,01 = \ x (t)\> \°i\ 2 and \ C z\ 2 are 
always equal to 1/2 and Ro is the center of the region 
where steps form. Moving away from this point, one of 
the |Cfc| 2 's becomes dominant (Fig. [3) and e gi (R,t) lies 
on top of the corresponding BOPES. 

To elaborate on how the TDPES switches between 
the two adiabatic states, we Taylor expand \C k (R, t)\ 2 
around Rq and keep only up to the linear order terms 



where 



«(*) = 



C l (R,t) =[l±a(t)(R-R )]/2, 



(V R \F 1 (R,t)\) Ra -^B.\F 2 (R,t)\) 1 
\X(t)\ 



(15) 



(16) 



Eq. jl5), using the relation < \C k \ 2 < 1 (k = 1,2) or 
equivalently \R — Rq\ < a" 1 , estimates the width of the 
region, AR, where the switching between BOPESs oc- 
curs with — , i.e., AR = —. Hence, the larger the values 
of a, the sharper the steps become. 

As an example, we discuss the TDPES at t = 31.87 fs 

in Fig. bl As it is seen, e gi switches from €~g {R) to 
(2) 

6 bo {R) over tri e region where \Fi \ and \F 2 \ cross (see the 
bottom plot). As I -Fx I and \F 2 \ have opposite slopes and 
cross where they are small, a is large yielding a small 
AR. Outside the switching region, one of the |Cfc| 2 s 
becomes dominant. Interestingly, the exact electron- 
nuclear density contains signatures of the behavior e gi , 

i.e., where e gi coincides with e% {R), presents one peak 

in analogy with |</>^(r)| 2 (see Fig. lj, while, where 

(2) 

e g i(R,t) follows t BO {R)t it displays two peaks, like 

l^i^( r )| 2 (Fig-[l} The step of e g i in the intermediate re- 
gion is indicated by scars in the full electron-nuclear 
density. 

In conclusion, we have presented generic features 
of the exact TDPES for situations in which, according 
to the standard BO expansion framework, significant 
non-adiabatic transitions occur and the nuclear wave- 
packet splits at the avoided crossing of two BOPESs. 
For the ID model system studied here, the TDPES is 
the only potential that governs the dynamics of the nu- 
clear wave-function (the vector potential can be gauged 
away) and provides us with an alternative way of vi- 
sualizing and interpreting the non-adiabatic processes. 
We have shown that the gauge-invariant part of the TD- 
PES, e g i(R,t), is characterized by two generic features: 
(i) in the vicinity of the avoided crossing, e gi (R, t), be- 




FIG. 3. Top: the full electron-nuclear density at the t = 
31.87 fs. Middle: a snapshot of the gauge invariant part of 
the TDPES (solid black line) at the t = 31.87 fs. For reference, 



,(!) 



■ rj0 (red dashed-line) and ej^ (green dashed-line) are shown. 
Bottom: Expansion coefficients (indicated in the figure) of the 
(two states) adiabatic expansion of the full wave-function and 
the exact electronic conditional wave-function (see the text) at 
the* = 31.87/s. 
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comes identical with a diabatic PES in the direction of 
the wave-packet motion, (ii) far from the avoided cross- 
ing, e g i(R,t), as a function of R, is piecewise identical 
with different BOPESs and exhibits nearly discontinu- 
ous steps in between. The latter feature holds after the 
wave-packet branches and leaves the avoided crossing. 
The gauge-dependent part, e g d(R, t), on the other hand, 
is piecewise constant in the region where e g i (R, t) co- 
incides with different BOPESs. Hence e g d(R-, t) has lit- 
tle effect on the gradient of the total TDPES, but may 
shift the BOPES-pieces of e g i(R, t) by different constants 
causing the exact TDPES to be piecewise parallel to the 
BOPESs. These features of the TDPES support the use of 
diabatic surfaces as the driving potential when a wave- 
packet approaches a region of strong NAC. Moreover, 
they are in agreement with the semi-classical picture 
of non-adiabatic nuclear dynamics that suggests to cal- 
culate the classical forces acting on the nuclei accord- 
ing to the gradient of only one of the BOPESs. We expect 
that these findings will ultimately lead to improved al- 
gorithms for the mixed quantum-classical treatment of 
electrons and nuclei. 
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